Polarimetric-interferometric boreal forest scattering model for BIOMASS end-to-end simulator

A polarimetric-interferometric forward model (FM) for extended covariance matrix modeling is presented. The FM has been designed to be used within the end-to-end simulator for BIOMASS, a new ESA satellite mission aiming at the global mapping of above-ground forest biomass with P-band synthetic aperture radar (SAR). The FM uses linear regression models for prediction of backscatter intensity and HH-VV correlation coefficient, and the random volume over ground (RVoG) model for the prediction of the interferometric correlation coefficients. For boreal forest, parameter values for these sub-models have been derived using polarimetric-interferometric SAR data acquired within the BioSAR 2007 campaign over the Swedish test site Remningstorp. The FM is evaluated qualitatively in a boreal forest scenario through a side-by-side comparison with BioSAR 2007 data. The general agreement is good, although there are regions with structures which cannot be reproduced by the model, probably due to insufficient forest description by the input parameters

Digital Canopy Model Estimation from TanDEM-X Interferometry using High-Resolution Lidar DEM

Interferometric TanDEM-X data are used together with high-resolution, airborne lidar-derived digital elevation models (DEMs) to produce digital canopy models (DCMs) for the boreal forests of Remningstorp and Krycklan, situated in southern and northern Sweden, respectively. An overview of interferometric data processing is given. First results showing the potential of TanDEM-X-based forest canopy mapping are presented. It is concluded that baselines giving height-of-ambiguity values in the order of 50-80 meters are preferable, although factors such as angle of incidence and along-track baseline are also of importance. Clear-cuts can easily be detected thanks to the high resolution of TanDEM-X imagery. Seasonal variations of scattering height are most visible for deciduous trees, where the scattering height is significantly lower in the winter, probably due to the lack of leaves. \ua9 2013 IEEE

Factorized Geometrical Autofocus for Synthetic Aperture Radar Processing

This paper describes a factorized geometrical autofocus (FGA) algorithm, specifically suitable for ultrawideband synthetic aperture radar. The strategy is integrated in a fast factorized back-projection chain and relies on varying track parameters step by step to obtain a sharp image; focus measures are provided by an object function (intensity correlation). The FGA algorithm has been successfully applied on synthetic and real (Coherent All RAdio BAnd System II) data sets, i.e., with false track parameters introduced prior to processing, to set up constrained problems involving one geometrical quantity. Resolution (3 dB in azimuth and slant range) and peak-to-sidelobe ratio measurements in FGA images are comparable with reference results (within a few percent and tenths of a decibel), demonstrating the capacity to compensate for residual space variant range cell migration. The FGA algorithm is finally also benchmarked (visually) against the phase gradient algorithm to emphasize the advantage of a geometrical autofocus approach

Regression-Based Retrieval of Boreal Forest Biomass in Sloping Terrain using P-band SAR Backscatter Intensity Data

A new biomass retrieval model for boreal forest using polarimetric P-band synthetic aperture radar (SAR) backscatter is presented. The model is based on two main SAR quantities: the HV backscatter and the HH/VV backscatter ratio. It also includes a topographic correction based on the ground slope. The model is developed from analysis of stand-wise data from two airborne P-band SAR campaigns: BioSAR 2007 (test site: Remningstorp, southern Sweden, biomass range: 10-287 tons/ha, slope range: 0-4 degrees) and BioSAR 2008 (test site: Krycklan, northern Sweden, biomass range: 8-257 tons/ha, slope range: 0-19 degrees). The new model is compared to five other models in a set of tests to evaluate its performance in different conditions. All models are first tested on data sets from Remningstorp with different moisture conditions, acquired during three periods in the spring of 2007. Thereafter, the models are tested in topographic terrain using SAR data acquired for different flight headings in Krycklan. The models are also evaluated across sites, i.e., training on one site followed by validation on the other site. Using the new model with parameters estimated on Krycklan data, biomass in Remningstorp is retrieved with RMSE of 40-59 tons/ha, or 22-33% of the mean biomass, which is lower compared to the other models. In the inverse scenario, the examined site is not well represented in the training data set, and the results are therefore not conclusive

Impact and modeling of topographic effects on P-band SAR backscatter from boreal forests

P-band SAR backscatter has been proven to be useful for forest biomass prediction. However, there is a need for further studies on effects of topography on P-band backscatter. In this paper, two prediction models for backscatter are evaluated, one using only biomass as predictor and one which also includes topographic corrections. Data from the BioSAR 2007 and BioSAR 2008 campaigns are used to evaluate the models. A multi-scale error model which is able to handle data from several imaging directions is used. For HH, the slope correction on stand level used in this paper is unable to correct for topographic effects. This is consistent with previous results that within stand topographic variability has a significant impact on HH P-band backscatter. For HV and VV, the model which considers topography gives lower prediction errors than the model which does not include topography. Moreover, for these polarizations topographic the correction strongly reduce the variability in backscatter measurements between imaging directions for stands with ground slopes larger than about 5 degrees

Topographic Correction for Biomass Retrieval from P-band SAR Data in Boreal Forests

The influence of the ground slope on radar backscatter has been proven to be greater for lower radar frequencies due to deeper canopy penetration. In this study, multiple heading, Pband SAR data of boreal forest in Sweden was used to find a model for topographic correction for improved biomass retrieval. Eleven models were tested and the best model was selected. The selected model was then used for biomass retrieval. Even by means of the most simplified approach, forest biomass could be established with a root-mean-square error of approximately 50 t/ha for HV and 66 t/ha for HH

Estimation of Forest Biomass From Two-Level Model Inversion of Single-Pass InSAR Data

A model for aboveground biomass estimation from single-pass interferometric synthetic aperture radar (InSAR) data is presented. Forest height and canopy density estimates Delta h and eta(0), respectively, obtained from two-level model (TLM) inversion, are used as biomass predictors. Eighteen bistatic VV-polarized TanDEM-X (TDM) acquisitions are used, made over two Swedish test sites in the summers of 2011, 2012, and 2013 (nominal incidence angle: 41 degrees; height-of-ambiguity: 32-63 m). Remningstorp features a hemiboreal forest in southern Sweden, with flat topography and where 32 circular plots have been sampled between 2010 and 2011 (area: 0.5 ha; biomass: 42-242 t/ha; height: 14-32 m). Krycklan features a boreal forest in northern Sweden, 720-km north-northeast from Remningstorp, with significant topography and where 31 stands have been sampled in 2008 (area: 2.4-26.3 ha; biomass: 23-183 t/ha; height: 7-21 m). A high-resolution digital terrain model has been used as ground reference during InSAR processing. For the aforementioned plots and stands and if the same acquisition is used for model training and validation, the new model explains 65%-89% of the observed variance, with root-mean-square error (RMSE) of 12%-19% (median: 15%). By fixing two of the three model parameters, accurate biomass estimation can also be done when different acquisitions or different test sites are used for model training and validation, with RMSE of 12%-56% (median: 17%). Compared with a simple scaling model computing biomass from the phase center elevation above ground, the proposed model shows significantly better performance in Remningstorp, as it accounts for the large canopy density variations caused by active management. In Krycklan, the two models show similar performance

Autofocus and analysis of geometrical errors within the framework of fast factorized back-projection

This paper describes a Fast Factorized Back-Projection (FFBP) formulation that includes a fully integrated autofocus algorithm, i.e. the Factorized Geometrical Autofocus (FGA) algorithm. The base-two factorization is executed in a horizontal plane, using a Merging (M) and a Range History Preserving (RHP) transform. Six parameters are adopted for each sub-aperture pair, i.e. to establish the geometry stage-by-stage via triangles in 3-dimensional space. If the parameters are derived from navigation data, the algorithm is used as a conventional processing chain. If the parameters on the other hand are varied from a certain factorization step and forward, the algorithm is used as a joint image formation and autofocus strategy. By regulating the geometry at multiple resolution levels, challenging defocusing effects, e.g. residual space-variant Range Cell Migration (RCM), can be corrected. The new formulation also serves another important purpose, i.e. as a parameter characterization scheme. By using the FGA algorithm and its inverse, relations between two arbitrary geometries can be studied, in consequence, this makes it feasible to analyze how errors in navigation data, and topography, affect image focus. The versatility of the factorization procedure is demonstrated successfully on simulated Synthetic Aperture Radar (SAR) data. This is achieved by introducing different GPS/IMU errors and Focus Target Plane (FTP) deviations prior to processing. The characterization scheme is then employed to evaluate the sensitivity, to determine at what step the autofocus function should be activated, and to decide the number of necessary parameters at each step. Resulting FGA images are also compared to a reference image (processed without errors and autofocus) and to a defocused image (processed without autofocus), i.e. to validate the novel approach further

Biomass Retrieval Algorithm Based on P-band BioSAR Experiments of Boreal Forest

A new biomass retrieval algorithm based on P-band multi-polarization backscatter has been developed and evaluated based on SAR and ground data over boreal forest. SAR data collections were conducted on three dates at a test site in southern Sweden (Remningstorp, biomass < 300 tons/ha; late winter to early summer 2007) and on a single date at a test site in northern Sweden (Krycklan, biomass < 200 tons/ha; fall 2008). The retrieval algorithm is a multiple linear regression model including the HV-polarized backscatter coefficient, the VV/HH backscatter ratio and the ground slope. Regression coefficients were determined from Krycklan data followed by algorithm evaluation using Remningstorp data. The results from the latter show that RMS errors vary in the range 29-42 tons/ha depending on date and stand type. The new algorithm is also compared with alternative algorithms and found to give significantly better performance. The developed model is a significant step towards an algorithm which gives consistent results across multiple sites and dates, i.e. when forest structure, topography and moisture conditions is expected to vary